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SUMMARY 


There  is  considerable  interest  in  the  subject  of 
efficient  wing  design  for  high  speed  flight,  together  with  good 
manoeuverability  at  transonic  speed  and  satisfactory 
near-to-the-ground  performance.  An  important  area  is  the 
detail  design  of  wing  leading  edges  with  or  without  leading 
devices  and  variable  camber,  the  flow  may  be  attached  or 
seoarated. 

In  this  report,  slender  bod''  theory  with  exact 
boundary  condition  is  used  to  calculate  the  flow  past  a  thin 
conically  cambered  wing  (i.e.  with  drooped  leading  edges). 

The  leading  edge  flow  separation  has  been  modelled  as  a 
vortex-cut  arrangement . 

Calculated  results  suggest  that  the  leading  edge  camber 
has  a  very  marked  effect  on  the  local  pressures  near  the  leaa- 
ing  edges,  and  hence  on  the  overall  aerodynamic  forces.  In 
particular,  the  lift-drag  ratio  can  be  enhanced  by  suitable 
choice  of  leading  edge  droop. 

An  interesting  feature  of  the  results  is  the  indication 
of  multiple-valued  solution  for  realistic  cambers,  i.e.  when 
the  shoulder  of  the  droops  is  verv  near  the  leading  edqe.  In 
experiment  this  is  not  likely  and  another  separation  svstem 
may  exist  near  the  high  curvature  wing  section.  The  work  in 
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this  report  therefore  points  to  inadequacy  of  current  design 
approaches  with  simulation  of  one  vortex  system  only. 

The  method  of  this  report  may  also  be  extended  to 
incorporate,  thin  extensions  of  wing  span,  leading  edge 
devices,  secondary  vortices  and  variation  of  separation 
points.  Cross-section  thickness  may  also  be  included.  More 
complex  vortex  sheet  representation  of  the  separation  may 
also  be  incorporated.  From  the  point  of  view  of  aircraft 
manoeuverability ,  wing  root  gaps  or  leading  edge  gaps  are 
of  interest.  Another  consideration  is  for  the  asymmetrical 
effects. 
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OA  in  Z-plane,  Fig. 5. 

radius  of  circle  in  ZQ-plane,  Fig. 5. 

aspect  ratio  =  4k. 

defined  by  equation  (33) . 

defined  by  equation  (32) . 

CB  in  Z-plane,  Fig. 5. 
c/a 

lift  coefficient, 
linear  lift  coefficient. 

non-linear  lift  coefficient. 

drag  coefficient. 

lift  induced  drag  coefficient. 

profile  drag. 

defined  by  Z  e^^  «  ej  +  i  fj 
function 

defined  in  Appendix  A. 
defined  in  Appendix  A. 
given  by  h2  =  1  +  (1/c2) 

defined  in  Fig. 4.  ~  amount  of  leading  edge  droop 

imaginary  part 
defined  by  eauation  (40) 
defined  in  Appendix  A. 
defined  in  Appendix  A. 

cot  (angle  of  sweep) 
lift  dependent  drag  factor  n  A  CD/CL2 
lift 
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M,N  defined  in  section  II. 7 

n  defined  in  Fig. 4. 

real  part. 

R  radius  of  circle  in  Z 3-plane  in  Fig. 5  (=  a/sin6cos6). 

r  OP  in  Zj -plane.  Fig. 5. 

r  CP'  in  Zi-plane,  Fig. 5. 

r'  parameter  in  equation  (10) 

s  semi-span 

T  vertical  slit  plane,  Fig. 5.  T  =  yT  +  iz? 

Te  vertical  slit  plane,  Fia.5.  Te  =  yxg  +  izT^ 

Tg  vortex  position  in  Tg-plane,  Fig. 5. 

l 

u,v,w  perturbation  velocities  in  x,y,z  system. 

V  velocity  at  infinity. 

vn  velocity  component  normal  to  surface  in  cross-flow  plane. 
W  complex  velocity  potential. 

Wj  (j  =  1,2,3, 4)  defined  in  section  II. 4. 

x,y,z  Cartesian  coordinates. 

Z  y  +  iz 

Zv  yv  +  iz^  vortex  position  in  Z-plane. 

Zj  yj  +  i  Zj  (j  =  0,1, 2, 3)  complex  planes  see  Fia.5. 

a  angle  of  incidence. 

aQ  angle  of  incidence  at  zero  lift, 

otp  defined  in  equation  (3?) 

a  angle  of  incidence  for  no  singularity  at  the  leading  edge, 

s 

B  tan  6 

6  half  the  angle  AO'B  in  Z2-plane,  Fig. 5. 

T  vortex  strength, 

e  droop  angle  as  defined  in  Fig. 4. 

defined  in  section  II. 6. 
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6  angle  POA  in  Z* -plane.  Fig. 5. 

k  lift  dependent  drag  factor  =  J1ACD/CL2. 

\  angle  PO'yo  in  ZQ-plane,  Fig. 5. 

p  density. 

<j>  half  angle  POB  in  Z2-plane,  Fig. 5. 

$  velocity  potential. 

•2-^  normal  component  of  velocitv  in  cross-flow  olane 

3n 

ip  defined  by  simp  =  tan4>/tan6 

suffix  u  refers  to  the  upper  surface. 

I  refers  to  the  lower  surface, 
superscript  '  refers  to  variable  quantity  in  integration 

subscript  L  or  L.E.  refers  to  leading  edge 
v  refers  to  vortex 

i 
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I  INTRODUCTION 


There  is  considerable  interest  in  the  subject  of  efficient  wing 
design  for  high  speed  flight  together  with  good  manoeuverabil ity  at  transonic 
speed  and  satisfactory  near-to-ground  performance.  An  important  area  is 
the  detail  design  of  wing  leading  edges  with  or  without 
leading  edge  devices  and  variable  camber.  The  flow  may  be 
attached  or  separated. 

The  designer  is  faced  with  simultaneous  constraints  of 
efficient  high-speed  cruise  (at  low  lift  coefficient  CL, 
high  Mach  number  M) ,  satisfactory  take-off  and  landing 
performance  (at  hicrh  C_  ,  low  M)  and  good  transonic 

ii 

manoeuverabil ity  (low  to  high  C L,  M  near  1).  These  constraints 
necessarily  result  in  performance  limiting  compromises.  By 
way  of  an  example,  we  take  the  case  of  an  aircraft  with 
highly  swept-back  wings  (Figure  la)  or  having  a  portion  which 
is  highly  swept-back  (Figure  lb) .  It  has  been  shown  that 
gains  in  flying  efficiency  i.e.  lift-drag  ratio,  can  be 
obtained  by  drooping  the  leading  edge.  The  amount  of  droop 
and  its  geometry,  e.g.  the  shoulder  position  depends  on  the 
operating  CL.  Thus  for  low  CL  values  the  optimum  position 
of  the  shoulder  of  the  droops  may  be  near  the  leading  edge, 
whilst  at  higher  values  of  C^,  better  results  are  obtained 
with  somewhat  gentler  camber  shapes  (Figure  lc).  The  choice 
of  droop  is  subject  to  compromises  of  this  nature. 
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The  extent  of  such  limiting  compromises  can  be  reduced 
by  permitting  the  designer  with  greater  flexibility  in  the 
geometric  configuration.  One  approach  is  the  use  of  fixed 
or  variable  geometry  near  the  leading  edges.  Variable  camber 
with  or  without  slats  or  other  devices  (Figure  Id)  fall  into 
this  category. 

For  conventional  aircraft  with  wings  of  low  sweep-bi 
there  is  a  considerable  amount  of  information  available  the 
subject  of  leading  edge  design  with  or  without  devices  (R  '  j.). 
Calculation  methods  have  been  devised  for  treating  2-D 
geometries  (Figs.  2  and  3).  The  methods  are  strictly 
applicable  only  when  3-D  effects  such  as  those  due  to  wing- 
tips  or  fuselage  junctions,  are  small. 

For  wings  of  higher  sweep-back,  however,  there  seems 
to  be  only  a  small  amount  of  information  available,  particularly 
with  regard  to  leading  edge  devices.  An  idea  of  the  order  of 
gains  from  leading  edge  devices  can  be  obtained  by  reference 
to  the  work  of  Ray  and  Hollingsworth  (Ref. 2)  on  F-4  Fighter 
aircraft  with  leading  edge  sweep-back  51.4°.  They  conclude 
that  incorporation  of  devices  resulted  in  a  sizable  33% 
improvement  on  the  buffet  onset,  I/D  performance  gain  of  35% 
at  CL  ®  .8  and  improved  lateral  directional  characteristics 
throughout  the  test  Mach  number  range  of  0.6  to  0.94.  The 
improvements  were  verified  in  a  subsequent  flight  evaluation. 

Goodmanson  and  Gratzer  (Ref. 3)  show  24%  improvement  in  L/D 
at  Cl  ■  0.4,  using  droop  and  slats  on  a  highly  sweep-back  wing. 
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At  present  there  are  no  suitable  methods  which  take 
into  account  the  flow  separation  and  are  generally  applicable. 
There  are,  however,  a  few  attempts  to  solve  the  3-D  wings 
with  separation  under  development.  These  are  of  panel  or 
collocation  type.  They  require  a  great  deal  of  computer 
time  and  use  large  cores  in  the  computer.  It  is  not  clear 
whether  these  methods  will  permit  adequate  resolution  with 
panels  or  collocation  points  to  deal  with  high  curvature  droops 
near  the  leading  edges,  or  leading  edge  devices,  or 
secondary  vortices. 

There  have  been  however  a  number  of  attempts  (Refs.  4  to 
16)  to  calculate  the  flowfields  within  the  framework  of 
slender  body  theory.  The  slender  body  theory  is  applicable 
in  subsonic  or  moderately  supersonic  flow  (i.e.  the  component 
of  flow  normal  to  the  wing  leading  edge  is  less  Mach  0.6). 

Refs.  4-7  concentrate  on  attached  flow.  In  particular, 

Cooke  (Ref. 7)  deals  with  very  general  and  realistic  droop 
configurations.  The  flows  calculated  are  however  for  the 
particular  case  where  the  leading  edge  singularity  vanishes. 
Refs.  8-12  idealise  the  flow  separations  by  line  vortex-cut 
arrangements.  Attempt  has  been  made  only  to  calculate  the 
flows  for  flat  or  circular  arc  camber  wings.  The  approach 
of  Jobe (10)  although  nominally  for  general  camber  ianores  the 
zero  incidence  flows. 

In  Refs.  13-16,  the  leading  edge  separations  are  more 
realistically  modelled  by  vortex  sheet  representations.  Once 
again  only  the  uncambered  wing  and  circular  arc  camber  cases 
have  been  fully  considered.  A  more  general  approach  by 


Fernandez  and  Holla  (16)  does  not  include  the  camber  effects 
at  zero  incidence  directly.  In  this  report  an  attempt  has 
been  made  to  calculate  the  flows  around  general  wing  camber 
shapes  with  droop  on  basis  of  Ref. 7  and  idealised  vortex-cut 
arrangements  as  in  Ref. 8.  This  procedure  is  considered  to 
be  an  essential  step  before  embarking  on  more  complicated 
vortex  sheet  representation. 

It  must  be  mentioned  that  the  attached  flow  solutions 
can  be  obtained  by  not  including  the  separation  effects. 

The  method  of  this  report  may  also  be  extended  to 
incorporate,  thin  extensions  of  wino  span,  leading  edge 
devices,  secondary  vortices  and  variation  of  separation 
points.  Cross-section  thickness  may  also  be  included. 


I 

I 
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II  THEORETICAL  FORMULATION 

II. 1  Equation  of  Motion 

Equation  of  motion  to  be  satisfied  is  the  well 
known  Laplace  Equation  which  represents  a  slightly 
perturbed  main  stream  velocity  V  and  corresponding 
Mach  Number  M : 


(1  -  M2)  $  + 


xx 


$  +$ 
yy  zz 


o. 


where  ♦  is  the  disturbance  velocity  potential  and 
x,  y,  z  are  Cartesian  co-ordinates  fixed  to  the  wing 
(x-axis  along  the  centre  line  of  the  wing  Fig. 4  ). 


The  flow-field  is  restricted  to  highly  swept  wings 
and  the  term  (1  -  M2)1?*^  may  be  neglected.  The 
equation  of  motion  then  becomes  Laplace's  equation  in 
y  and  z  dimensions 


4>  + 

yy 


i>  =  o. 

zz 


This  is  the  governing  equation  of  motion  of  slender 
wing  (or  body)  theory  and  applies  throughout  the  Mach 
number  range. 


II. 2  Boundary  Conditions  and  Flow  Model 
The  conditions  on  the  wing  are: 

(i)  The  wing  is  solid  and  hence  the  normal  velocities 


are  zero. 


(ii)  For  attached  flow  at  the  leading  edge  the  condition 
is  that  the  flow  turns  around  the  leading  edge.  If 
the  edge  is  sharp  the  velocity  is  infinite. 

If  the  flow  separates  at  the  leading  edge,  the 
condition  to  be  applied  is  that  the  velocity  is  finite  at 
the  edge. 

The  conditions  in  the  field  are: 

(i)  The  disturbances  vanish  at  infinity. 

(ii)  The  fluid  Pressure  is  continuous. 

A  conical  wing  with  leading  edge  drooped,  at  incidence  o  to  the 
freestream  V  and  with  Cartesian  axes  xyz  is  shown  in 
Fig. 1.  The  semi-span  of  the  wing  is  s=kx  and  the 
shoulder  is  distance  ns  from  the  centre,  the  amount  of 
droop  is  Hs. 

The  model  approximates  the  vorticity  in  the  separated 
vortex  sheets  from  the  leading  edge  (±  s(x) ,  -Hs)  by  a 
concentrated  pair  of  vortices  of  strength  r (x)  above  the 
wing  at  positions  ±yv<x),  zy(x).  The  strength  r <x) 

varies  in  streamwise  direction  and  therefore  in  order  to 
satisfy  Kelvin’s  Law  of  Conservation  of  Circulation,  feeding 
vortex  sheets  in  the  form  of  'cuts'  of  strength  ^ ^ 
at  the  leading  edges  have  been  assumed.  The  axis  of 
vorticity  in  these  sheets  is  assumed  to  lie  in  the  yr-plane 
and  because  slender  wing  conical  flow  is  implied,  the  cuts 
do  not  affect  the  velocity  profile  in  cross-flow  plane. 
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It  is  noted  that  the  component  of  free  stream  V 
in  x-direction  is  Vcosct  and  in  the  cross-flow  yz-plane, 
Vsinct.  Both  components  must  be  considered  to  obtain 
the  velocity  and  pressure  fields. 


II. 3  Camber  Shape  Transformations  and  their  Properties 

Ir.  order  to  make  a  satisfactory  study  by  the  slender 
body  theory,  it  is  necessary  to  find  conformal  trans¬ 
formation  which  transforms  the  camber  line  in  the  cross- 
flow  plane  Z  =  y  +  iz  into  a  circle  or  a  vertical  slit. 
For  a  general  shape  this  cannot  be  done  in  closed  form 
and  in  the  present  approach,  the  series  of  transformations 
as  used  by  Maskell  and  Cooke- Ref  .7  are  employed.  These 
are  summarised  as  follows  (see  also  Fig. 5): 


Transformations 
Z  =  y  +  iz 

Z2  =  z2  —  4c2  or 

Z 2  ~  Z s  +  2iatan{  or 

Zi  -  f2  t  \h\  -  4a 2  or 

Z  =  Z 1  -  iatan 6  or 

o  1 


T  ■  T  +  2ia.sin6  or 
e  » 


Z2  =  +  4c2  (1) 

Z  3  =  Z  2  -  2iatan6  (2) 

Z,  -  Z,  +  (3) 

Z,  =  Z  +  iatanS  (4) 

o 


zo  *  r  1  +  4a>:  <5> 

Te  -  T  -  2  ia^sinti 


(6) 


These  transformations  leave  the  point  at  infinity 
unchanged . 

It  must  be  mentioned  that  the  camber  in  Z-plane 
is  generated  on  the  assumption  of  a  circular  arc  in  the 
Z 3 -plane. 

The  Transformation  derivatives  are  given  by  the 
following  relations 


We  shall  also  require  for  the  evaluation 

of  Boundary  Condition  of  zero  force  on  Vortex-cut 
arrangement  (Section  II. 6)  and  this  is  derived  here. 


Prom  equations  (7), 


we  have 


dZ  dZ  dZ2  dZQ 

dT  dZ3  *  dZi  ‘  dTe  * 


d2Z  /dZ  d2Z  dZ 3  d2 Z  dZ  d2Z0  dT 

*  ’  dTr/dT  "  dZ 3dT  *  dZ  dZidT  *  dZ2  dT7  *  dZ0 


1  dZ ,  1  dZ  2a2  dZ!  dZi 
Z3  *  dT  Z  '  dT  Z]!'  dT  dZ2 


+ 


2a|  Z 
(ZQ+ai 


T 


(8) 


If  a  point  P  on  the  circle  is  given  by 
ZQ  =  a  sec6ei?l  (=  aie1*)  and  Zj  =  rei6' 
then  it  can  be  easily  shown  that 

r2  -  a2  =  2r.a. tanfisinO  (9) 

r  =  a  tan3sin6  +  a/1  +  tan26sinz0  (10) 

and  that 

y3  =  2a  cos0  /I  +  tan26sin*0 
z3  -  -2a  tan6cos20. 

If  we  write 

sin<!>  *  sin6cos0  (11) 


then  we  have 
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where 

Hence 


yj  *  R  sin  2<j> 
z 3  =  -  R(l-cos24>) 

R  =  2a  cosec  26. 

Z3  =  2R  sin<)>e-i^ , 


(12) 


and  0  is  angle  shown  in  Fig. 2  and  R  is  the  radius  of 
the  circle  in  the  Z 3 -plane. 


Family  of  Cross  Section  Shapes 
There  are  simple  relations  between  the  various 
parameters  of  the  drooped  wing  and  these  are  derived  here. 

At  the  leading  edge  of  the  wing  A  where  0=0 
and  <J>  =  6,  Z3  and  Z  -  (+s  -  iHs)  are  related  by 

_  j  r 

Zj  s  2a  sec6e  . 

Z  =  (s  -  iHs)2  =  4c2  +  4a2sec26e  2i<S. 

On  separating  real  and  imaginary  parts,  we  get 

s2  (1  -H1 )  =  4c2  +  4a2  (1  -  tan26)  (13) 

s*H  =  4a2tan6.  (14) 

At  the  shoulder  point  B  where  Z3  =  0  and 
y  *  2c  =  ns.  (15) 

From  equations  (13),  (14)  and  (15)  we  find 

1  -  H2  -  n2  1-tan2  5 
H  "  tan* 

or 

cot  26  =  1  -  ~ 

cH 


Given  n  and  H  this  determines  6  and  then  c/a 


may  be  obtained  from 

-i  c2  n2tan<S 

c  *  I2 - —  * 

It  is  useful  at  this  stage  to  collect  here  for 
reference  certain  formulae  which  will  be  required  later. 

If  ds^  is  an  element  of  the  arc  of  the  circle,  then 
o 

dso  =  a  secfidX  =  {r2  +  (||)2}*d0  =  £^d0*  (16) 

For  points  on  droops  we  have  from  equations  (1)  and 

(12) 

2  Z  (|J  +  i||)  =  8R2  sin<)>e* 3  i<J> 
or 

dy  idz  _  4R2sln4> 

d$  d<t>  "  * 


If  we  let 

Z  e3i<^  =  ei  +  ifi 
then  on  the  droops 

dz  f  i 
dy  =  "  eT 


(17) 


Now  if  ds  is  the  element  of  arc  on  the  droops  we  note  that 


ds  =  dy  ,  .dz 
c34>  d$  d$ 


4R2  sln4> 


(18) 


on  the  droop 


I 

1 

I 

I 


Since  the  surface  of  the  wing  is  conical  through  origin 
0  its  equation  must  be  homogenous  in  x,y,z  and  may  be  written 
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h  =  f(fe> 

and 

|f  -  k  f(y/kx)  -  £  f*  (y/kx) 

^  vz _ ky  .  dz 

kx  kx  dy 

=  fs,(eiz  +  fiz) 

using  equation  (17)  and  putting  kx  =  s. 

Hence 

If  =  fii  { (e5  +  if>)(y  +  iz>> 


The  velocity  normal  to  the  contour  in  the  cross-flow 
plane  due  to  the  component  V  cosoi  is  on  the  "droops" 
(ref. 7),  see  Fig. 6. 


il 

3n 


V  cosaf^  V  cosa.^-  .  (Z2e3i^) 

v  =  -  _2*_  =  - -  sel  - 

n  {1  +  (f§)2}*  /I  +  (f,/ei)2 


_  k  V  cosa  (Z2e3i<^) 

•  •  i  •  • 

•1*1 

vn  is  zero  on  the  flat  part  of  the  section. 


(23) 


In  order  to  find  the  complex  velocity  potential  we  shall 
need  the  component  v^  in  ZQ-plane  which  is  related  to 
vn  by  the  mapping  ratio  j^|  . 

(19),  (21)  and  (23)  we  find  that 


i 


Using  equations  (1) ,  (12) 
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16k  V  a2cosa  (c2sin3<f>  +  R2sin34>)  sln6cos8 . 
s  cos2 6  I Z | 2r 


. .. (24) 


Care  must  be  taken  with  signs  here.  We  adopt  the  convention 

that  vR  represents  the  component  ot  velocity  along  outward 

drawn  normal.  Therefore  vn  is  positive  on  upper  surface 

and  negative  on  the  lower  surface  of  the  droops.  vnQ  is  zero  at 
the  leading  edge. 


II. 4  The  Complex  Velocity  Potential 


The  Complex  Velocity  potential  W  comprises  the  cross- 


flow  and  axial  flow  contributions.  The  cross  flow 


contributions  arise  due  to  the  freestream  components  V  sina, 
V  cosa  and  the  vortices  ±T  at  (±y^,  z^) •  The  axial 
contribution  arises  due  to  V  cosa.  This  can  be  represented 


as 

W  -  W|  (V  sina)  +  W2  (V  cosa)  +  W3  (D  +  W*  (x) .  (25) 

By  virtue  of  transformation  of  the  Z-plane  to  Te-plane, 

W|  and  Wj  become 


Wi 

Wj 


-i  V  sina.Te 


ir 
"  Iff 


log 


(26) 


(27) 


-  22 


where  Te^  is  the  position  of  right  hand  vortex  in  the 
Te-plane. 


Wt  is  simply  written  as 


W*  =  V  cosa.x. 


(28) 


W2  arises  for  wings  other  than  thin  and  flat  and  Weber  (Ref. 17) 
and  Cooke  (Ref. 7)  show  that  it  is  best  represented  in 
lQ-plane.  W2  is  defined  as  the  complex  velocity  potential 
required  to  produce  a  normal  velocity  vnQ  at  the  surface 
of  the  circle.  It  is  equivalent  to  considering  the  effects 
of  a  source  of  strength  211 .  aj  .  at  the  centre  of  the 
circle  together  with  a  source  distribution  on  the  circumference 
of  strength  2^^  -  v^)  per  unit  length,  where  is 

the  mean  velocity  on  the  circumference  given  by 


v'  (A')dX' 

O 


(29) 


For  a  thin  wing  it  can  be  verified  that  v^  vanishes  as  expected. 
Hence  the  complex  velocity  potential  W2  is  given  by 


*2 


2v^  log  (Z0-Z')ds' 


a.  f2* 

Jo  <*'>  log  (Zo-Z')  d*'  (30) 


where  ZQ'  is  a  point  on  the  contour  but  Z0  is  any  point  in  the 
complex  plane. 
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On  the  surface  of  the  circle  in  ZD -plane,  following  Weber^17\  the 
velocity  potential  *2  can  be  written  as 
a  f2ir 

♦2  =  (W2)  =  -jp  I  Vn  (X')  log  | Zq-Zq I  dX'  (31) 

*  ft 


we  shall  also  need 


dW2 

dZ^ 


for  calculation  of  velocities.  It  follows  from 


differentiation  of  equation  (30) 
•2ir  v ^  (X")  dX" 

o 


f2  aj  f2 

dz° '  *  L 


Zo  -  Zo 


*v  <ze> 


(32) 


The  reduction  of  integral  Iv  (ZQ)  is  dealt  with  in  Appendix  A.  This 
is  interpreted  as  a  Cauchy  Principal  Value  integral  and  its  value  depends 
on  whether  the  point  Zq  lies  on  the  surface  of  the  circle  or  outside  it, 
as  follows 

iX 


Zq  on  circle  ZQ  =  aje 

r2ir  v 


dW2 


_  x  ^  V  _  ai  f2ir  V"o  (X'}  '  vnQ  (X)  , 
~  Iv  <Zo>  -  J  Z~ ^  d 


X'  + 


iX 


Zq  outside  circle  Z0  =  pe  ,  p  >  a  1 


ft.  1. 


dZ, 


(ZQ) 


-  f 
*  L 


2ir  v  '  (X") 


vn_  (X) 


a,  vn  (X) 
1  llo 


2a,  v_  (X) 


(33) 


Zo  -  Zo 


dX'  + 


(34) 


dW, 


Alternative  forms  for  W2  and  its  derivative  have  been  derived  by 
,  (7)  dZo 

Cooke  .  He  uses  9'  as  the  integration  variable  instead  of  X'.  Using 

equations  (16),  (24)  and  (30),  W  can  be  written  as 

r2ir 


r* 

W2  =  C  B'  (0')  log  (Z0  -  ZT0)  d9' 


where 

16  k.V  a2  cos  o 


C  = 


ITS  cos2  6 


B'  (0')  =  +  R2sin3$*)  sinB^cosO* 

cos$"| Z'|z 

*  On  the  surface  of  the  wing  W2  can  be  reduced  to  #2  and  evaluated 
exactly  (Ref.  7  and  17)  as  shown  in  Appendix  B. 

dW2 


(35) 

(36) 

(37) 


The  complex  velocity  derivative  is  given  by 


dW2  _  „  f2w  B'  (8')  d 
‘  Jo  Zo  -  % 


d0' 


Equation  (20)  on  substitution  of  equations  (21),  (22)  and  (23)  becomes 

i£ 

2ir 


X  -T 

W  (Te)  =  -  iv  sina  .  Te  -  ~  log  ~  ♦  W2  +  V  cosa.x 


with  W2  given  by  either  equation  (30)  or  (35), 


(38) 


( 19) 
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dW 

The  complex  velocity  derivative  is  given  by 


dW 

dZ 


dW 

dT- 


dW 


dTe  _  ...  .  dTe  ir  f  1 

'  dZ  lVsino  .  dz  ~  2n  (.Te-Tei 


Te+Tei 


dTe  dW2  d^o 
"dZ"  +  dZ^  *  "dZ" 


(40) 


with  -r=—  given  by  either  equation  (33)  or  (34)  or  (38). 
a^o 

II .5  The  Boundary  Condition  of  Finite  Velocity  at  the  Leading  Edge 

This  condition  is  applied  in  the  cross-flow  plane  at  the  leading  edge 

— i6  dW 

(Z  =  s  -  i  Hs,  ZQ  =  a3e  ).  The  complex  velocity  derivative  •£=  is  non 

singular  at  the  leading  edge.  Using  equation  (39),  we  obtain 

dZ. 


dW  _  .  dTe  ir  f  1 

dZ  "  lVsino  •  dZ  '  2ir  iTe-Tei 


Te+Te J 


dT.  dW, 
~dZ~  +  dZ^ 


o 

dZ 


(41) 


where 
dZ„  dZ 


“o 

dZ 


_ o 

dZj 


dZj 

dZj 


dZ„ 


dZ3  dZj 


dZ3  *  dZ  dZ2 


dz3 

dZ 


^e 

dZ 


dTe 

dZ, 


dZ, 


o 

dZ 


At  the  leading  edge  d2l  is  infinite,  although  other  derivatives  of  the 


must  have 


r  dw ' 


dZ2 

fini 

=  0 


_  dVJ 

transformations  are  finite.  Hence  for  finite  ^  at  leadinS  edg®»  we 


L.E. 


f-£-)  =  -£[t4 - M  ’  5?  ’  3rl  =  ° 

( dTeJ  l.E.  L  2  (Te  T®1  Te+Te^  dZ°  dTdL.E. 


or 


.  „  .  if 
l  V  sina  = 


Te  j  +Te i 


m—  +  “D 


h  ei 


fdZQ] 

^Jl.E. 


(42) 


where 
aD  - 

and 


dW- 


dZr 


L.E. 


dZo) 

Zr2 

•dTe^L.E. 

L.E. 


Op  can  be  obtained  directly  from  equation  (33)  by  numerical  integration 
or  from  the  exact  method  outlined  in  Appendix  C. 
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II. 6  The  Boundary  Condition  of  Zero  Force  on  the 
Vortex-cut. 

The  condition  of  zero  force  on  the  vortex-cut  combination 

follows  directly  from  Brown  and  Michael  (8)  and  can  be  written 

as  v 

(V,  +  iw*)Z=Zv  =  x  (2ZV  ‘  V 

or 

(V,  -  iw,)z=2v  -  V  <«> 

where  vj  +  iwx  is  the  velocity  at  the  right-hand  vortex  and 
is  found  by  subtracting  the  velocity  field  of  the  vortex  at 
Zv  from  the  total  velocity  and  taking  the  limit  as  Z  -*■  Zy. 

Using  equation  (38)  ,  the  complex  velocity  potential  Wj 
at  the  right-hand  vortex  is  written  as 

Wi  =  W  +  log(Z-Zv)  (4  4) 


and  by  differentiating  with  respect  to  Z  and  taking  the  limit 
as  Z  -*■  Zv,  we  obtain 


i  •  e  • 


iVsinct  -  |£(-i- 
2u  T-T 
e  6! 


iT  1 
211  Z-Zy  * 


dW2 


with  ^r-  given  by  either  equation  (34)  or  (37). 
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Substituting  equation  (35)  for  V  sina,  we  have, 

.  irVei  _  dZ^,.  .  iT  fj  1  §T»  1  1  T 

vt-iwi»  “2nT"^Te1  D5^ "  2n  [iTg-T^" 

/dW2  dZp\ 

+  \dZo  ’  dzA_ 


The  complex  equation  (46)  now  gives  a  complex  value  for 
the  vortex  strength  for  a  given  position  of  the  vortex 
The  procedure  we  follow  is  to  fix  the  height  of  the  vortex  and 
vary  the  spanwise  position  until  r  has  zero  imaginary  part. 
The  angle  of  incidence  then  follows  directly  from  equation  (42). 
This  procedure  is  based  on  Reference  8. 
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II. 7  Pressure  Distribution  and  Forces 
Pressure  Coefficient 

The  pressure  coefficient  to  the  same  order  as  the  linearized  equation 
for  the  potential  is  given  for  slender  configurations  by 

Cp  =  sin2a  -  2  ~  cosa  -  i-r-  ($  2  +  ♦  2  ) 

v  v2  y  z 

where 

*  =  *  (x,  y,  z)  JQ«. 

In  conical  flow  ♦  may  be  written  as 

*  =  x.#  +  y.$  +  z.$  +  U  (a  constant) 

x  J  y  z 

The  differentiation  with  respect  to  x  is  for  constant  y  and  z,  i.e.  constant 
Z,  so  that 

V*  -  %j '  *  9  [f j- 

It  follows  then  that 

Cp  =  sin2a  -  f[[^(W)  -  -  0[SjJ«s  a 

■  7  [RS)J  •[>«]•]•« 

dH 

W  and  ^  follow  from  equation  (39)  and  (40).  The  constant  v  must 
be  chosen  so  that  Cp  vanishes  as  y,  z  -*■  <*>. 


On  the  surface  of  the  wing  the  part  of  the  complex  potential 


due  to  drooped  leading  edges  (W2)  presents  some  problem  and  care 
must  be  taken  in  its  evaluation.  In  Appendix  B,  a  method  for  reduction 
of  W2  to  #2  based  on  Refs  7  and  17  has  been  given. 


Similarly,  care  must  be  exercised  in  evaluation  of  .  The  Cauchy 
Principal  value  integral  is  implied  and  numerical  integration  procedure 


should  account  for  this. 
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The  lift  force  coefficient  based  on  the  projected 

area  of  the  win?  s2/k  is  given  by  combining  the  linear  lift 
expression  of  Ref. 7  with  the  non-linear  lift  due  to  the 
vortices  (e.g.  Ref. 15) 

C 

—  *  nk2  .  4  (§)  (2+2c2+tan2<$)^££ 

cos 2 a  S  k 

-  4<f>’§§f£j  •  +  2  •  TC  •  (Tei%,>< 

where  M  -  1  +  8c2  +  16c"  -  16|2  -  , 

N  =  2  (5c2  +  85“  -  40c6  -  14|2  +  12|"  +  40^'  +  |s). 


The  lift  and  induced  drag  coefficients  and  may  also 

be  obtained  from  integration  of  pressures  on  the  wing  as  follows 
Forces 

We  define  the  conical  winr,  surface  by  the  equation 


F(x,y,z)=  z  -  kx  f  (^/^x)  =  0 

Unit  normal  vector  £  defining  the  surface  is  given  by 
i  F  +  i  F  ♦  k  F 

a  =  —2 — — =-£ 


where 


<r„2  *  ryz  ♦  r72)-' 


=  -  k  f  <y/kx)  .»/,!•  <y/kK) 

*  -  k  e/kx>  .  k  (y/kk)  Jk/dy 


r,  *  -  k,  f'  <y''k*>'^  =  <'  <y'kx>  *  d%y 


F*  *  1 


We  also  define 


F  =  (F  2  +  F  2  +  F  z}-‘ 
r  x  y  z 
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1 


The  Force  vector  dF 
dF  =  Cp 

where  d  S  is  an  elemental  surface  area.  This  is  can  be  represented  in 
terms  of  arc  length  d  o  as 

tJS  s  i  ,  x  .do  for  a  conical  wing  of  centre  line  chord  x. 


is  defined  by 
d  S  n 


The 


do  =  dy2+dz2  =  dy2  (1  ♦  g)  =  dy2  U  ♦  f'2> 
force  coefficients  C  and  C.  in  x  and  z  directions  respectively  are 


then  given  by 


y=: 


i  dF  2 


y=-r 


y=s 


Cp.dfl.  *x 
F„ 


y=-r. 


am) 


.  [ cp  1-»>V/,C)  d  ,y/kx, 

■'-]  r 

.  y=s  y=r. 

2_  (  >  k  dr  ?  (  Cp  df,  £z 

=  17  »v=-:.  '  Kx5  y=-::  'V 


2  l  *  k  dr  .  ? 

kx'  ’  y=- 

cv  wry-  „  ,v/fcx, 


Lift  coefficient  and  drag  coefficient 

coefficient*  C  and  C  . 

X  z 

C.  =  C  cos  a  -  C  sin  « 

u  Z  X 

Cn  -  C  sin  a  +  C  cos  a 

U  Z  X 


follow  by  resolving  the  force 


It  must  be  mentioned  that  although  the  above  derivations  are  generally 
applicable,  in  this  report  wr  have  not  included  viscous  effects. 
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III  RESULTS 


The  type  of  camber  lines  generated  by  the  method  are 
shown  in  Fig.  7.  The  angle  of  droop  measured  at  the  leading- 
edge  increases  as  the  shoulder  position  moves  outboard. 

For  each  configuration,  there  are  two  particular  incidences 

art  and  a  .  a  corresponds  to  zero  lift  and  a  to  the 
o  s  o  s 

attitude  where  the  load  at  the  leading  edge  vanishes.  For 

an  uncambered  wing  a  and  a  both  coincide  at  0°.  However, 

for  a  wing  with  leading  edge  droop,  both  a  and  a  will 

0  s 

be  positive. 

Strictly  speaking  the  value  of  aQ  depends  on  the  method 
of  solution  of  the  problem  i.e.  if  the  flow  is  assumed  to  be 
separated  as  in  present  theory  or  if  it  is  attached.  The 
difference  in  practice  for  "realistic"  configurations  is 

however  likely  to  be  very  small.  a  therefore  may  be 

o 

calculated  easily. 

The  angle  a  has  a  special  significance,  it  indicates 
the  side  of  the  wing  on  which  the  leading  edge  vortex  lies. 

For  a  >  ag  the  leading  edge  vortex  lies  above  the  upper 
surface  and  vice-versa. 

Fig.  8  depicts  the  variation  of  sinog/k  for  a  set  of 
camber  lines.  As  might  be  expected,  aQ  increases  both  with 
increasing  amount  of  droop  and  as  the  shoulder  position  approaches 
the  leading  edge. 


As  mentioned  in  the  Introduction,  previous  work 
contains  adequate  reference  to  uncambered  and  circular  camber 
wings  (zero  shoulder  position) .  The  emphasis  in  this  report 
is  therefore  on  the  aerodynamic  effects  as  the  shoulder 
position  moves  out  toward  the  leading  edge.  Both  attached 
and  separated  flow  solutions  have  been  considered. 
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III.l  Attached  Flow 

In  general ,  the  Attached  Flow  theory  applied  for  small 
incidences  only.  Fig.  9  shows  the  effect  of  shoulder 
position  on  at  fixed  incidences  for  two  values  of  camber 

parameter  H  ■  0.05s  and  0.1s.  It  is  noted  that,  in  general, 
as  incidence  increases  the  effect  due  to  a  fixed  droop  decreases. 
When  the  shoulder  of  the  droop  tends  to  the  leading  edge,  there 
is  a  gain  in  lift,  although  for  most  part  droop  causes  a 
reduction  in  lift  at  a  particular  attitude. 

In  Fig.  lo  pressure  distributions  are  shown  at  5° 
incidence  (k  =  .25)  for  various  shoulder  positions.  It  is 
interesting  to  note  the  oscillatory  behaviour  of  the  pressures 
near  the  leading  edge  on  the  upper  surface  as  the  shoulder  moves 
outboard  from  ,775s  to  .875s. 

Such  a  behaviour  is  also  likely  to  arise  for  some  other 
droop  configuration  at  a  different  incidence.  It  also 
illustrates  the  problem  of  determining  CL  and  from  Cp 

distribution  by  numerical  integration.  A  large  number  of 
sampling  points  may  be  required  near  the  leading  edge. 


:i 

ii 

•j 
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III. 2  Separated  Flow 

Figs.  11  and  12  show  a  selection  of  results  as  the 
shoulder  position  of  the  droop  is  varied.  The  camber 
parameter  values  considered  are  H  =  .05s,  .10s.  in  each 
case  the  following  relationships  have  been  presented: 

(i)  vortex  locus. 

(ii)  vo  *tex  height  (2V  +  H)/&  against  a. 

r 

(iii)  vorcex  strength  ^  against  o<  . 

(iv)  CL  against  a 

Constant  or  a  lines  where  appropriate  have  been 

interpolated. 

For  shoulder  positions  less  than  ,9s,  we  infer 

(i)  vortex  locus  moves  outwards  both  with  increasing 

shoulder  position  and  camber. 

(ii)  For  a  given  shoulder  position,  the  vortices  move 
closer  to  the  surface  as  camber  increases. 

(iii)  For  a  given  camber,  as  the  shoulder  moves  outwards 
the  vortex  strength  decreases  at  a  fixed  incidence. 
This  is  particularly  noted  at  smaller  incidences. 

At  higher  incidence  the  curves  approach  the  flat 
wing  case. 

(iv)  For  a  given  camber,  as  the  shoulder  moves  outward, 
the  lift  decreases  at  a  particular  incidence,  the 
tendency  being  more  marked  at  lower  incidences. 
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An  interesting  feature  of  results  is  that  for  shoulder 
positions  very  near  the  leading  edge,  the  vortex  height  is 
apparently  multi-valued  through  a  small  incidence  range 
Ce.g.  for  configuration  with  n  «  .95s,  H  =  .05s  and 
10.0°  <  a  <  15.0°).  Lift  '''incidence  and  vortex  strength  ~ 
incidence  relationships  also  exhibit  a  similar  tendency. 

The  solution  therefore  admits  the  possibility  of  up  to  three 
different  flow-fields.  Fig. 13  showsthe  three  types  of 
pressure  distributions  corresponding  to  three  vortex  positions  - 
"low",  "intermediate"  and  "high"  at  a  £  11.7°.  The  span- 
wise  variation  of  the  velocity  tangential  to  the  upper  surface 
in  the  cross-flow  plane  is  shown  in  Fig. 14.  for  all  three 
vortex  positions.  The  curves  indicate  that  reattachment 
streamlines  which  enclose  the  separation  and  then  split  in 
two  ways  impinge  on  the  wing  for  the  "low"  and  "intermediate" 
vortex  positions  (vis.  .98s  and  .85s  respectively)  but  not 
for  the  "high"  vortex  position.  Flow  patterns  as  sketched 
in  Fig.  15  are  indicated.  The  flow  patterns  also  show 
general  agreement  with  the  features  of  pressure  distributions. 

For  the  "low"  vortex  position  the  peak  suction  appears  outboard 
of  the  shoulder.  For  the  "intermediate"  and  "high"  vortex 
position  the  peak  suctions  occur  near  and  inboard  of  the  shoulder 
position  respectively. 

This  analysis  leads  to  the  question  of  what  happens  in 
experimental  flows  where  the  results  are  not  likely  to  be 
multiple-valued  for  a  given  configuration  and  there  would  be  a 
preferred  flow  field.  Two  possible  explanations  may  be  advanced 


as  follows:- 


-  35 


(i)  The  preferred  flow-field  would  depend  on  the  shape  of  the 

wing  plan  form  near  the  section  under  consideration.  If  the  sweepback 
decreases  aft  of  the  section  under  consideration  (as  in  Fig. 16  (a)) 
the  tendency  may  be  towards  attached  flow  type  of  flow  field 
e.g.  either  for  "low"  or  "intermediate"  position.  Conversely, 
if  the  sweepback  increases  aft  of  the  section  under  consideration, 
then  the  tendency  may  be  for  "high"  position  flow  field. (Fig.l6(b)) 

(ii)  The  velocities  on  the  wing  surface  for  the  "low"  and 
"intermediate"  vortex  positions  are  fairly  large  near  the 
leading  edge.  It  is  therefore  not  difficult  to  visualise  the 
existence  of  multiple  vortex  system  as  sketched  in  Fig.  17. 

Of  the  above  two  explanations,  the  latter  one  seems  very 
plausible  and  experimental  evidence  indicates  the  presence  of 
a  number  of  vortex  systems.  This  however  suggests  that  the 
theoretical  approaches  generally  used  for  design  work  which 
simulate  only  one  separation  system  are  not  likely  to  be 
adequate.  This  criticism  applied  equally  to  conical  and 
non-conical  or  lifting  surface  approaches. 

It  must  be  mentioned  that  Levinsky  and  Wei  (17)  also 
show  the  existence  of  multiple  separations  on  slender  bodies 
with  Strakes  (Fig. 18).  Their  model  however  deals  with  only 
one  separation.  Additional  criteria  based  on  the  development 
of  the  boundary  layer  on  the  body  will  be  required  to  fix  the 
strength  and  position  of  the  second  separation. 
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Comparisons  at  Constant  C^. 

It  is  of  interest  to  compare  various  quantities  for  a 
given  CL.  Taking  for  example  a  value  of  CL  =  0.3,  we 
look  at  the  effect  of  droop  shoulder  position  keeping  the 
height  of  camber  a  constant  at  H  =  .05s.  The  results  are 
shown  in  Fig.  19  and  20  we  note  the  following  as  droop 
shoulder  position  increases: 

(i)  linear  part  if  lift  increases  and  the  non-linear 
part  decreases. 

(ii)  induced  drag  decreases  and  then  increases  again. 

(iii)  incidence  required  at  the  centre  line  increases. 

(iv)  The  vortices  move  closer  towards  the  leading  edge. 

(v)  Up  to  shoulder  position  of  0.85s,  the  peak  suction 
pressure  on  the  upper  surface  increases  and  also  its 
position  moves  outwards  signifying  a  drag  reduction 
acting  on  the  drooped  leading  edge. 

For  shoulder  position  0.95s,  the  suction  peak  is 
much  smaller  and  this  corresponds  to  a  slight  increase 
in  lift  induced  drag. 

To  enable  general  conclusions  to  be  drawn,  this  type  of 
analysis  needs  to  be  carried  out  for  a  number  of  droop 
configurations. 


Lift  Dependent  Drag  Factor  kj 


it  A  Cd£ 

The  variation  of  lift  dependent  drag  factor  kj  =  — —  g—  a-t  various  lift 

'-'L 

coefficients  for  a  representative  series  of  wings  with  leading  edge  droop  has 
been  considered. 

Wings  with  Circular-arc  Cambers  n  =  0 

Figure  21  shows  the  variation  of  the  factor  kj  for  leading  edge  droop  up 
to  0.4s.  It  is  noted  that  the  factor  reduces  with  increasing  Cl.  The 
discontinuities  in  the  curves  near  low  Cl  values  occur  near  the  attachment 
angles  of  attack.  The  values  for  kj  without  leading  edge  singularity  based 
on  Smith's  results  are  shown  for  reference.  Flow  separation,  therefore, 
is  beneficial  and  gives  a  reduction  in  kj. 

(14) 

Figure  22  shows  a  comparison  between  Barsby's  results  with  vortex 
sheet  separation  model  and  the  present  approach.  It  is  interesting  to  note 
that  at  lower  values  of  CL,  the  present  approach  gives  lower  values  for  kj 
but  for  higher  values  it  predicts  higher  values.  The  correspondence  of  the 
two  methods  improves  with  increasing  camber. 

Wings  with  Leading  Edge  Droop 

The  variation  of  shoulder  position  for  leading  edge  droop  of  0.05s  and 
0.1s  has  been  depicted  in  Figures  23  and  24.  It  is  noted  that  moving  the 
shoulder  of  the  camber  outwards  means  a  reduction  in  kj. 

For  a  given  shoulder  position  at  0.8s  the  effect  of  leading  edge  droop 
has  been  shown  in  Figure  25.  It  is  noted  that  increased  droop  is  generally 
beneficial,  but  it  may  lead  to  a  limiting  value  for  droop.  Further  work 
can  be  done  on  this  aspect  to  optimise  the  leading  edge  droop  geometry  for 
given  Cl» 


IV _ CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  report,  slender  body  theory  with  exact  boundary 
condition  has  been  used  to  calculate  the  flows  past  a  thin 
conically  cambered  wing  (i.e.  with  drooped  leading  edges) . 

The  flow  separation  has  been  modelled  as  a  vortex-cut 
arrangement . 

Calculated  results  suggest  that  leading  edge  droop  has 
a  very  marked  effect  on  the  local  pressures  near  the  leading 
edges  and  hence  the  overall  aerodynamic  forces.  In 
particular,  the  lift-drag  ratio  can  be  enhanced  by  suitable 
choice  of  leading  edge  droop.  The  calculations  need  to  be 
generalised  for  a  set  of  cambers. 

The  results  indicate  a  multiple  valued  solution  for 
realistic  cambers  i.e.  when  the  shoulder  of  the  droop  is  very 
near  the  leading  edge.  In  experiment,  this  is  not  likely  and 
another  separation  system  may  exist  near  the  high  curvature 
wing-section.  The  work  in  this  report  therefore  points  to 
inadequacy  of  current  design  approaches  with  simulation  of  one 
vortex  system  only. 

The  method  of  this  report  may  also  be  extended  (see  Fig. 
26)  to  incorporate  thin  extensions  of  wing  span,  leading  edge 
devices,  secondary  vortices  and  variation  of  separation  points, 
cross-section  thickness  may  also  be  included.  More  realistic 
vortex  sheet  representation  (Ref. 13)  may  also  be  incorporated. 

From  the  point  of  view  of  aircraft  manoeuverability ,  wing 


i 


root  gaps  or  leading  edge  gaps  may  also  need  to  be  studied. 
Asymmetrical  configurations  are  also  of  interest. 
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APPENDIX  A 


Evaluation  of  Integral  I  (Z  ) 

_ _ _  v  o 


From  Equation  (27)  we  have 
I  <2  I  .  ll  f  '"o  <»'>  d»' 

*  »  *  L  — r=!- — 


o  o 


,2* 

a,  I  V"«(A',"vn  <A>  a.  f2w  ,  , 

=  •  l  W  dA'  "  T  v"°(*>  {  4=5^ 


-o 

.  al  f  "  vn0(A',‘V"o(A> 

*  ~  JL  - r^T - dA  +  a»  Vno(A>  •  UZJ 


2n 


U- 

O  o 


u  o 


It  is  noted  that  the  value  of  I  (Z  )  depends  on  whether  Z  lies  on 

u  o  o 

the  surface  of  the  circle  or  outside.  It  can  be  evaluated  as  follows 
2* 
o 

2it 

a  ,  a 

r-z  ~x  dz 

o 

dA'  i 


dA 

z^T 


W  *  7  j 

O 

.  1  (‘’  1  .  dA  .  r  . 

-  •  J  d t*  c,'o 


o  o  o 
iA 


How  Zq  -  a^e  therefore  wr  nee  ^ 


7. 


o 


W  =  -  7  C^Kr>:': ' 5  ‘  ^  /‘Tr*  -  rV]dV 

o  o  n  o  /.  I_o  vo  O  J 


If  Z_  lies  on  the  snrfac-e  o f  the  circle  then 
o 

1  i Z  )  =  -  (2*i  -»i)ri 

U  O  w  /#  t, 

o  e 


If  Z  lies  outnide  the  circle  tln'ii 
o 

W-  --f  <2"-°>  -  l 

o  o 
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APPENDIX  B 


The  reduction  of  W?  to  $2  and  evaluation. 


(211 

W2  =  C  B'  log  (Z  -Z^dQ'. 
J0 


Denoting  the  real  part  of  W2  by  *2  we  have 

f2n 

*2  =  C  B'  log  |z0-z£|de*. 
j0 


In  the  ZQ-plane  if  $2  is  to  be  evaluated  at  a  point 
outside  the  circle,  the  integral  presents  no  problem. 
However,  on  the  surface,  special  care  is  required  to  deal 
with  log  IZq-Z^I  term. 


If  we  denote  the  radius  of  the  circle  by  d  then  we 

have 

♦  ,  -  c  j^B  log|y0+lz0~y0'-iz0'  i  lYo+iz0+^o'~izo'  ^ 

-J h 

from  the  symmetry  about  the  zQ-axis.  Hence  on  the  circle 


♦2  ■ 


C 

|  J  P  log 
-IJr2 


(yo-y0,)2  +  (zo-zo’>2 


(y0+yo')2  +  (*o-*o,) 


c 

1 


J^B'log  4|(d2-z0z0')2  -  y02y0,2}de' 
-ifc 


Putting  yQ2  *  d2-z02  we  have 


i 
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f1^2  .  i 

*2  *  c  B‘ log  2d | z0-z0' |d0 ' 


on  the  surface  of  the  section. 


Since  z0-z0'  =  zi-zi'  we  have  on  the  surface,  taking 
note  of  the  signs  of  B' 


*2 


t^/2 

C  B' log 
J0 


Z  1-Z  ' 
z  l“Zl 


hi 

I 

Jl 


d0' 


(Al) 


where  Z|u  is  value  of  z^  at  a  point  specified  by  0' 
and  zx  '  is  its  value  at  the  corresponding  point  on  the 

X/ 

lower  surface  specified  by  -6'. 


Given  sin\J>  = 


tant 

tan<5 


sin  4>  =  sin6cos0 


we  have 


COS0  = 

cos  $  = 

G2  = 


sinij/ 

G 

cos6 

G 


sin0  = 

sin$>  = 


1  -  cos2ipsin2  <5 . 


cosi|;sin6 

G 

sinipsimS 

G 
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r  —  a  tan6sin0  a/1  +  tan*  6s:* 

=  a/G.(l  +  sin6cosi|>) 

r  =  a/ g .  ( 1  "  sin6cosi(0 
2 

rr  =  a * 
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we  have 


and 


=  r  sin© 


a  cos6cosi h 
l-cos^sin6 


r'  sin  fl* 


a  cos 6 cost];1 
l-cos^'sin6 


z\  =  -r 'sin© ' = 
At 


a  cos6cosi P' 
1+cost^ '  sin<$ 


d© 


COS3  COS  Ip 

cos6sin© 
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On  substitution,  the  log  term  in  equation  (  Al)  becomes 


log 


zi-2}u 
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log 
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(1 -cosies  in  6)  (1-cosijr  sln6) 
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*2  now  becomes 
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where 
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If  B'de'  is  now  put  in  the  form 


B'd0 '  = 


-  jtan6cos26  g  * )  d^;  * 


where 


then 
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u 

-  jtan6cos26  J  g(^’)  log 
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Now  g (II~xp * )  =  -g(i|)') 


therefore 
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Ifi  may  also  be  reduced  in  similar  way 
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l+cos^ 1 sin6 
1-cos’jj '  sin  <5 
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$2  is  now  given  by 


*2  -  C 


j|-tan6cos26  (JA(ip)  +  JB) 
fn 
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The  value  of  JA(i{0  depends  on  it>  hut  is  a  constant 

for  a  given  wing  and  may  be  calculated  by  numerical  integration. 
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Evaluation  of  J,  (i|<) 
fn 

*  j  g  (tp ' )  log  |  cos’f'-cosH'  *  j  d\J> 1 . 


This  evaluation  of  thin  form  of  integrals  has  been 
discussed  by  Weber  (17)  and  a  solution  using  Fourier  Series  is 
suggested.  We  know  that  g(4>)  is  an  odd  function  which  also 
gives  g(<Ji'=0)  *  g(\J>'=Ifc)  =  g(ijj'=n)  =  o  we  write 

N-l 

g(\l>')  =  I  bvcosv4»'  (N  even) 

v=l ,3,5 


N  j>H 
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N  K 

=  ji  y  — ^  cos  v  ^ . 

V»1  v 


The  coefficients  bv  are  calculated  from 
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At  the  leading  edge  ZQ  *  •  zo"zo  =  Zj-Z{  =»  a  - 
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Where  r'  is  value  of  r'  when  8'  is  replaced  by 
we  find  that 


Using  equations  (9),  (10)  and  (11),  and  noting  that 
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Using  equation  (36)  we  obtain 


a  =  3  2Vi  cosa.e1^.— a*  XF.  j 
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Equation  (35)  becomes 
T  +T  . 
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The  integral  I,  has  been  evaluated  in  Ref. 7  as 
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